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What is Claimed is : 

1. A method of estimating the spatial variation of 
magnetic resonance imaging radio frequency (RF) signal 
intensities within an object from measured RF signal 
intensities of a uniform spin density medium surrounding the 

5 object, said method comprising: 

acquiring a magnetic resonance image of an object 
bounded by a medium which is of essentially uniform spin 
density, on the length scale of resolution of the image; 

formulating a semi-empirical mathematical model of the 
10 spatial variation in RF signal intensity within the said 
object; and 

fitting said model to selected measured RF signal 
intensities of said surrounding medium to obtain an estimate 
of the spatial variation of magnetic resonance imaging 
15 radiof requency (RF) signal intensities within said object. 

2. A method according to claim 1 wherein said step of 
formulating said semi-empirical mathematical model of the 
spatial variation in RF signal intensity within the object 
comprises : 

5 locating a plurality of points in a plane of said image 

which are notionally considered to act as apparent receivers 
of RF signals (hereinafter referred to as "RF receiver 
points' 7 ) ; and 

formulating a spatial profile of RF signal intensity 
10 within said object in said plane of said image relative to 
the said RF receiver points. 

3. A method according to claim 2 wherein said RF receiver 
points are located on said surrounding medium. 

4. A method according to claim 2 wherein said spatial 
intensity profile is formulated in a manner which provides 
a concentric reduction in RF signal intensity with 
increasing distance from said RF receiver points. 



-25 - 



5. A method according to claim 4 wherein said reduction 
in RF signal intensity is formulated as having the same 
rate of reduction in signal intensity with increasing 
distance from all said RF receiver points. 



6. A method according to claim 5 wherein said spatial 
profile is formulated as having an exponential reduction in 
signal intensity with increasing distance from said RF 
receiver points. 



7. A method according to claim 5 wherein said spatial 
profile is formulated as having a reduction in signal 
intensity with the reciprocal of the distance from said RF 
receiver points raised to a chosen power. 



8 . A method according to claim 6 wherein 

when the magnetic resonance image is acquired in an axial 
plane said semi-empirical mathematical model is formulated 
as : 

T/ . T -Rj(x-x„) 2 [\/f.cos 2 0 n +fsm 2 0 n \(y-y„) 2 \l/f.sm 2 0„+fcos 2 0 n 'l-(x-x n Xy-y„)[\/f.sin 20„ -/sin2 <9„ ] 

where x and y are image coordinates, I(x,y) is the model of 
the spatial intensity profile of the RF signal, N is the 
number of RF receiver points, I n is the estimated signal 
intensity at the n th RF receiver point and is to be 

10 determined by the fitting procedure, R is the rate of decay 
in signal intensity with distance from the RF receiver 
points and is to be determined by the fitting procedure, x n 
is the x position of the n th RF receiver point, y n is the y 
position of the n th RF receiver point, f is the degree of 

15 ellipticity to the concentric reduction in RF signal 
intensity with increasing distance from the RF receiver 
points and is to be determined by the fitting procedure, and 
0 n is the angular alignment of the major axis of the 
elliptical contour to the decay in the RF signal intensity 
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20 from the n tb RF receiver point with the coronal plane. 

9. A method according to claim 6 wherein 

said semi-empirical mathematical model is formulated as: 

where x and y are image coordinates, I(x,y) is the model of 
5 the spatial intensity profile of the RF signal, N is the 
number of RF receiver points, I n is the estimated signal 
intensity at the n th RF receiver point and is to be 
determined by the fitting procedure, R is the rate of decay 
in signal intensity with distance from said RF receiver 

10 points and is to be determined by the fitting procedure, x n 
is the x position of the n th RF receiver point, y n is the y 
position of the n th RF receiver point, and f is the degree 
of ellipticity to the concentric decay in RF signal 
intensity with increasing distance from said RF receiver 

15 points and is to be determined by the fitting procedure. 

10. A method according to claim 9 wherein said selected 
measured RF signal intensities comprise pixels of maximal 
signal intensity in respective radial line segments of an 
image of the surrounding medium (hereinafter referred to as 

5 "selected pixels") . 

11. A method according to claim 10 further comprising 
connecting said selected pixels by a line to form a trace. 

12. A method according to claim 11 further comprising said 
method further comprises deriving from said trace local 
maxima in signal intensity for each of said radial line 
segments . 

13. A method according to claim 12 wherein when said 
magnetic resonance image is acquired in an axial plane said 
RF receiver points are located at the local maxima in signal 
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intensity of the said selected pixels. 

14. A method according to claim 13 wherein when said 
magnetic resonance image is acquired in said axial plane the 
said local maxima in signal intensity of the selected pixels 
are no closer than 12 pixels. 

15. A method according to claim 14 wherein said local 
maxima in signal intensity determined from the selected 
pixels serve as starting values for the signal intensities 
I n of the said RF receiver points in the said fitting 

5 procedure. 

16. A method according to claim 15 wherein said fitting 
procedure includes deriving an initial estimate of the rate 
R of decay of signal intensity with distance from said RF 
receiver points determined by: 

5 locating local minima amongst said selected pixels in a 

manner consistent with determination of said local maxima; 

for each RF receiver point locating two nearest local 
minima and fitting an exponential signal intensity decay 
curve as a function of the distance from the said RF 

10 receiver point to the signal intensities of the two local 
minima constrained to the signal intensity at the position 
of the said RF receiver point to determine an estimate of 
the rate R n of decay in signal intensity with distance away 
from the said RF receiver point; and 

15 calculating the average of the estimates of the rate of 

decay R n in signal intensity with distance from each RF 
receiver point which is to be used as an initial estimate 
of R. 

17. A method according to claim 9 wherein when the 
magnetic resonance image is acquired in the axial plane the 
elliptical contours of the decay in the RF signal intensity 
from a given RF receiver point have their minor axes 

5 approximately normal to a plane of an RF receiver coil 
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element of an MRI machine acquiring said image of said 
object in closest proximity to the given RF receiver point. 

18. A method according to claim 11 wherein when the 
magnetic resonance image is acquired in the axial plane the 
angular alignment of each of the elliptical contours of the 
decay in the RF signal intensity from the RF receiver points 
5 is determined relative to said trace. 

19.. A method according to claim 18 wherein the angular 
alignment of the major axes of the elliptical contours of 
the decay in the RF signal intensity from the RF receiver 
points is determined relative to said trace by: 
5 centering a window kernel on each RF receiver point and 

finding the two furthest positions from said RF receiver 
point at which said trace of the selected pixels passes 
through the perimeter of the window kernel; and 

calculating an angle of a line bisecting said two 
10 furthest positions relative to the coronal plane. 

20. A method according to claim 19 wherein said window 
kernel which is centered on said RF receiver points to 
determine the angular alignment of the elliptical contours 
of the decay in the RF signal intensity from said RF 

5 receiver points is a square window kernel not less than 13 
pixels in width. 

21. A method according to claim 9 wherein said model 
spatial intensity profile I(x,y) is fitted to the image 
intensities of the selected pixels for optimum 
parameterisation of I n , R and f to obtain an estimate of the 

5 spatial variation in RF signal intensity throughout the 
object . 
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22. A method according to claim 9 wherein said method 
further comprises determining a RF spatial attenuation 
profile with which to rescale the magnetic resonance image 
intensities within the said object and thus minimise the 

5 spatial variation in RF signal intensity throughout the 
object . 

23. A method according to claim 22 wherein said RF spatial 
attenuation profile is determined by dividing the estimated 
RF signal intensity profile I(x,y) by the minimum or maximum 
value of I(x,y) within the region occupied by the object. 

24. A method according to claim 23 wherein the spatial 
variation in RF signal intensity throughout the object is be 
minimised by dividing the image intensities within the 
object by the RF spatial attenuation profile. 

25. A method according to claim 1 further comprising when 
the magnetic resonance image is a spin echo image , providing 
an estimate of the image intensities within the object at 
zero echo time by dividing the estimated signal intensities 

5 within the object at the given echo time by a percentage 
decay in signal intensity expected within the surrounding 
medium at the given echo time and then dividing this result 
by a ratio of hydrogen proton spin density expected of the 
bounding medium relative to hydrogen proton spin density 
10 expected of the object itself. 

26. A method according to claim 25 wherein calculating 
transverse relaxation rates within the object from a series 
of spin echo images acquired of the object and surrounding 
medium at different echo times. 

27. A method according to claim 26 wherein when the 
magnetic resonance image is a spin echo image said method 
further comprises using the estimates of the image 
intensities within the object at zero echo time in the 
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5 calculation of transverse relaxation rates within the object 
where the estimates are determined from a spin echo time 
image of shortest duration in the spin echo image series. 

28. A method according to claim 1 wherein when the object 
under consideration is the liver said surrounding medium is 
a layer of subcutaneous fat surrounding the abdomen. 

29. A method according to claim 28 wherein when transverse 
relaxation rates are calculated within the liver from a 
series of spin echo images the subcutaneous fat is used to 
estimate image intensities within the liver at zero echo 

5 time to assist in the calculation of the transverse 
relaxation rates. 

30. A method of enhancing a magnetic resonance image (MRI) 
of an object comprising: 

obtaining an initial MRI image of the object; 

estimating the spatial variation of magnetic resonance 
5 imaging radio frequency (RF) signal intensities within the 
object from measured RF signal intensities of a uniform spin 
density medium surrounding the object by: 

acquiring a magnetic resonance image of an object 
bounded by a medium which is of essentially uniform spin 
10 density, on the length scale of resolution of the image; 

formulating a semi-empirical mathematical model I(x,y) 
of the spatial variation in RF signal intensity within the 
said object; 

fitting said model to selected measured RF signal 
15 intensities of said surrounding medium to obtain an estimate 
of the spatial variation of magnetic resonance imaging 
radiof requency (RF) signal intensities within said object; 
and 
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rescaling the initial MRI image by dividing the signal 
20 intensities of the image of the object. 



